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Abstract 

We have used a tight-binding model that is fit to first-principles electronic-structure calculations 
for titanium to calculate quasi-harmonic phonons and the Gibbs free energy of the hexagonal 
close-packed (hep) and omega (u) crystal structures. We show that the true zero-temperature 
ground-state is the to structure, although this has never been observed experimentally at normal 
pressure, and that it is the entropy from the thermal population of phonon states which stabilizes 
the hep structure at room temperature. We present the first completely theoretical prediction of 
the temperature- and pressure-dependence of the hep-w phase transformation and show that it is 
in good agreement with experiment. The quasi-harmonic approximation fails to adequately treat 
the bec phase because the zero-temperature phonons of this structure are not all stable. 

PACS numbers: 63.20.Dj, 64.30.+t, 64.70.Kb, 65.40.Gr 



The experimentally observed ground state of titanium at room temperature is the hexago- 
nal close-packed (hep) structure |l| , even though zero-temperature first-principles electronic- 
structure calculations j3, 0, 0, Ojls | predict the omega (u) structure to be the energetically 
favored ground state. There are two possible causes for this discrepancy: (1) this is an- 
other failure of the local-density approximation (LDA) for calculating electronic-structure 
of solids, or (2) thermal effects stabilize the hep phase. In this paper we show that the second 
argument is correct: the thermal occupation of phonon states favors the hep structure above 
a pressure-dependent transition temperature. Only by including the vibrational entropy in 
the first-principles calculations can theory be brought into agreement with experiment. 

Titanium is a good case for testing our ability to predict solid-solid phase transformations 
with first-principles methods. It is an element that displays a rich phase diagram with several 
recently discovered high-pressure . When alloyed with other elements, it leads 

to materials of major technical importance for the airline, space, and other industries. 

It is possible to use first-principles methods to calculate the zero-temperature internal 
energy, &q(V), as well as the free energy contributions from the ions, Fi(T,V), and the 
electrons, Fe(T, V), to give the complete equation of state 

F(T, V) = %(V) + Fj(T, V) + F E (T, V). (1) 

Because of the large computational effort required to calculate the free energy of Eq. [T]over a 
significant range of volumes and temperatures, we have chosen to use a tight-binding model 
instead of conventional first-principles electronic-structure methods. We will demonstrate 
the high accuracy of our model and believe that it should be essentially identical with the 
best available first-principles methods. 

The model's accurate emulation of first-principles calculations (in this case the Wien97 



code 
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]) requires optimization of a single tight-binding parameterization such that the model 
closely reproduces the relative total energies of selected crystal structures. The structures are 
chosen to be both simple and representative: simple enough to make the first-principles cal- 
culations feasible, and representative of all microscopic environments the model is aimed at 
describing. A comparison of total energies as a function of volume with structures fulfilling 
these requirements are shown in Fig. Not shown but included in the fit as indepen- 
dent structures are fee crystals distorted to correspond to the longitudinal and transverse 
phonons at the reciprocal-space high-symmetry point X. These additional structures fine- 



tune the model's ability to correctly describe lattice vibrations [lo|. Crucial for convergence 
of the fitting procedure is to fit the cubic-structure energy bands at high-symmetry points 



in reciprocal space 
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The electronic contribution to the free energy depends on the electronic density of states 
(DOS) as a function of volume, n(E, V). The occupation of these states, given by the Fermi 
distribution f(E,T) = [e(- E ~ E /)/( fcsT ) + determines their entropy[13], 

S d (T,V) = -k B J[f\nf + (l-f)\n(l-f)}n(E,V)dE, (2) 

and hence the electronic contribution to the free energy F E (T, V) = —TS e i(T,V). This 
contribution is small compared to that of the ions but is included for completeness. 

The ionic contribution to the free energy depends on the phonon DOS, g(u, V), through 
the zero-point energy, 

U zero {V) = \ [ hwg(u,V)du, (3) 
as well as the temperature-dependent, quasi-harmonic free energy, [13J] 

F H (T, V) = k B T [ ln[l - e-^ kBT }g(to, V)du. (4) 
Jn 

This differs from the volume-independent harmonic approximation in that the quasi- 
harmonic approximation uses phonons calculated for small displacements (i.e., harmonic) at 
each volume. The resulting phonon frequencies are therefore volume dependent but temper- 
ature independent (they are calculated at zero temperature). The approximation requires 
that the phonons for any given crystal structure and volume are stable at zero temperature, 
which may not be true if the crystal can lower its total energy by spontaneously and continu- 
ously distorting into another crystal structure. At elevated temperatures (near melting) this 
approximation also neglects some important anharmonic contributions to the free energy 
that can be responsible for stabilizing some high-temperature phases. 

We have used the tight-binding model to calculate the electronic and the phonon DOS by 
evaluating the relevant energy eigenvalues and dynamical matrices on a fine mesh of wave 
vectors in the first Brillouin zone. In both cases the spectrum is smeared with a Gaussian, 
chosen to be as small as possible while keeping the DOS smooth and continuous. The mesh 
is refined until the free energy converges. 

The dynamical matrix at a given wave vector q is the Fourier transform of the force 
constants, which we calculate from the tight-binding model by the direct-force method [nl 
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171 ]. This method requires large simulation cells consisting of repeated unit cells 



transposed by vectors £, the forces on all atoms are calculated in response to the displacement 

s contain 54 (hep), 81 (u>), and 128 atoms 



id]. 



of the atom(s) in one unit cell. Our simulation eel 
(bec); for more details we refer the reader to Ref. 

Figure El compares our calculated phonon dispersion for the hep structure to experimental 
data 18]. The overall agreement is good, in particular near the zone center (around T) and 
along the vertical direction (parallel to the c-axis). Frequencies of modes in horizontal direc- 
tions (perpendicular to the c-axis) away from the zone center are somewhat high compared 
to experiment. The main characteristics, which the tight-binding model captures, give a 
reliable DOS, the relevant entity used in Eqs. 01 and |U to calculate the phonon contribution 
to the free energy. 



The phonon density of states is calculated for the hep 



231 ] . u, and bec structures at a 



dozen volumes corresponding to pressures in the range from zero to 170 GPa and is then 
used to evaluate the Gibbs free energy Q, G(T, P) = F(T, V) + PV . 

Comparison of the Gibbs free energies as a function of pressure and temperature for 
the hep and u crystal structures allows us to identify the phase transformation line. This 
approach is extremely demanding because tiny errors in the relative free energies can dramat- 
ically alter the phase boundary. Figure El shows good agreement between our theoretically 
calculated a — > uj phase transformation line with the experimental phase diagram^|. The 
calculated transition temperature at zero pressure is approximately 280 K. 

At normal pressure we know of no experimental data that shows that the u phase is 
the correct zero-temperature ground state. However, it has recently been shown that energy 
barriers exist for this phase transformation [20] that may make it extremely difficult to verify 
at low temperatures. Thus the kinetics may be just too sluggish for the system to bring 
itself into true thermodynamic equilibrium. 

We were not able to similarly calculate and include in Fig. El the transition into the 
bec structure. The open bec structure appears to be dynamically stabilized by entropy 
and anharmonic effects. At zero temperature some phonon modes are unstable due to a 
spontaneous instability of the bec crystal structure into the u phase. There is no energy 
barrier to this continuous transformation. In particular, this causes the longitudinal phonon 
mode with wave vector q = |(1, 1, 1), which distorts the bec crystal into the u phase, to 
be unstable below pressures of roughly 40 GPa. The mechanical instability of this phonon 



mode at low temperatures is also seen experimentally; the martensitic phase transition 
from the bcc phase to the hep phase prevents the high-temperature bec phase from being 
quenched to room temperature j^]. One branch of transverse phonon modes with wave 
vectors (£, £, 0), which are related to the bcc to hep transition, also appears unstable in our 
treatment. Both of these instabilities are consistent with strictly first-principles calculations 
of zero-temperature phonons 22]. 

Figure E^i shows the relevant segment of the (£,£,£) branch. The energy of the corre- 
sponding pathway from the bcc phase to the u phase, shown in Fig. develops a local 
minima for the bcc phase at pressures of roughly 40 GPa, but at lower pressures the crystal 
can deform from bcc to u without encountering an energy barrier. This pathway can be 
viewed as planes moving in accordance with the wave vector q = | (1,1,1) or as chains 
of atoms along the body diagonal moving with respect to each other. Calculations of the 
energy needed to move such a single chain of atoms indicate that for small displacements 
the movement is energetically favored, whereas the movement of a single atom is not. 

If only the stable phonons are included in the evaluation of Eq. 01 the bcc structure is 
calculated to be lowest in Gibbs free energy at pressures that are an order of magnitude 
higher than experimentally observed. A correct calculation requires a more sophisticated 
treatment that includes phonon-phonon interactions, i.e., anharmonic contributions to the 
free energy. 

In conclusion, we have for the first time used a completely theoretical approach, without 
any experimental input, to calculate the Gibbs free energies of competing crystal phases 
in titanium. Comparison of the hep and uj Gibbs free energies results in a pressure- 
and temperature-dependence of the structural phase transition in good agreement with 
experiment and serves to explain how finite-temperature entropy effects stabilize the ex- 
perimentally found room-temperature hep crystal structure over the theoretically found 
zero-temperature w-phase ground state. 

We thank Matthias Graf, Carl Greef, and Dallas Trinkle for helpful and encouraging 
discussions. This research is supported by the Department of Energy under contract W- 
7405-ENG-36. All FLAPW calculations were performed using the Wien97 package jjj. 



[1] S. K. Sikka, Y. K. Vohra, and R. Chidambaram, Prog. Mat. Sci. 27, 245 (1982). 



[2] M. J. Mehl and D. A. Papaconstantopoulos, Europhys. Lett. 60, 248 (2002). 
[3] J. S. Gyanchandani, S. C. Gupta, S. K. Sikka, and R. Chidambaram, J. Phys. Condens. Matter 
2, 301 (1990). 

[4] R. Ahuja, J. M. Wills, B. Johansson, and O. Eriksson, Phys. Rev. B 48, 16269 (1993). 
[5] G. Jomard, L. Magaud, and A. Pasturel, Philos. Mag. B 77, 67 (1998). 
[6] C. W. Greeff, D. R. Trinkle, and R. C. Albers, J. Appl. Phys. 90, 2221 (2001). 
[7] Y. K. Vohra and P. T. Spencer, Phys. Rev. Lett. 86, 3068 (2001). 
[8] Y. Akahama, H. Kawamura, and T. L. Bihan, Phys. Rev. Lett. 87, 5503 (2001). 
[9] P. Blaha, K. Schwarz, and J. Luitz, WIEN97, A Full Potential Linearized Augmented Plane 
Wave Package for Calculating Crystal Properties (Techn. Universitat, Wien, Austria, 1999). 
[10] S. P. Rudin, M. D. Jones, C. W. Greeff, and R. C. Albers, Phys. Rev. B 65, 235114 (2002). 
[11] D. A. Papaconstantopoulos, Handbook of the Band Structure of Elemental Solids (Plenum, 

New York, 1986). 
[12] M. D. Jones and R. C. Albers, Phys. Rev. B 66, 134105 (2002). 
[13] R. K. Pathria, Statistical Mechanics (Pergamon Press, Oxford, 1972). 
[14] K. Kunc and R. M. Martin, Phys. Rev. Lett. 48, 406 (1982). 
[15] S. Wei and M. Y. Chou, Phys. Rev. Lett. 69, 2799 (1992). 
[16] W. Frank, C. Elsasser, and M. Fahnle, Phys. Rev. Lett. 74, 1791 (1995). 
[17] K. Parlinski, Z. Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 78, 4063 (1997). 
[18] C. Stassis, D. Arch, B. N. Harmon, and N. Wakabayashi, Phys. Rev. B 19, 181 (1979). 
[19] D. A. Young, Phase Diagrams of the Elements (University of California Press, Berkeley, 1991). 
[20] D. R. Trinkle, R. G. Hennig, S. G. Srinivasan, M. D. J. D. M. Hatch, H. T. Stokes, R. C. 

Albers, and J. W. Wilkins, cond-mat/0301084. 
[21] W. Petry, A. Heiming, J. Trampenau, M. Alba, C. Herzig, H. R. Schober, and G. Vogl, Phys. 

Rev. B 43, 10933 (1991). 
[22] K. Pcrsson, M. Ekman, and V. Ozolins, Phys. Rev. B 61, 11221 (2000). 
[23] The c over a ratio of the hep structure is varied to give the lowest zero-temperature energy, 
i.e., in particular at high pressures the calculation is done for a hexagonal crystal with a c/a 
larger than the ideal value. 





10 14 18 10 14 18 

Volume (A 3 ) 

FIG. 1: Calculated energies for the crystal structures in the fitting database. Symbols are the 
first-principles results, solid lines are calculated from the fitted tight-binding model. Also included 
in the fit (but not shown here) are the fee crystal with distortions corresponding to the longitudinal 
and transverse phonons at the reciprocal-space high-symmetry point X. The uo structure lies lowest 
in energy, in contrast to the experimentally found ground state, hep. 
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FIG. 2: Experimental (crosses) [lg] and tight-binding (solid lines) phonon dispersion of titanium 
in the hep crystal structure at ambient pressure. Although some experimental details are not 
reproduced by our tight-binding model, the overall agreement is quite good. The phonon DOS 
that is relevant for the calculation of the Gibbs free energy should still be highly accurate. The 
square at T is from a first-principles frozen-phonon calculation. 




FIG. 3: The phase diagram of titanium. The dashed lines connect the experimental data points 
given by Young[l^|: the solid line shows our calculated a — > u> transformation. Thermal stabiliza- 
tion explains why at room temperature the experimentally observed hep structure is favored over 
the to structure (the calculated zero-temperature ground state). 




FIG. 4: (a) Calculated zero-temperature longitudinal phonon modes around the wave vector 
q = |(1,1,1), which distorts the bec crystal into the u phase, (b) Calculated zero-temperature 
energy landscape for the bec to u> pathway. 



